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Starting from a simple microscopic model, we have derived a kinetic equation for the matrix 
distribution function. We employed this equation to calculate the conductance G in a mesoscopic 
F'/F/F' structure with a domain wall (DW). In the limit of a small exchange energy J and an 
abrupt DW, the conductance of the structure is equal to G^d = ^o~-\o~i/{o~-\ + o{)L. Assuming 
that the scattering times for electrons with up and down spins are close to each other we show 
that the account for a finite width of the DW leads to an increase in this conductance. We have 
also calculated the spatial distribution of the electric field in the F wire. In the opposite limit of 
large J (adiabatic variation of the magnetization in the DW) the conductance coincides in the main 
approximation with the conductance of a single domain structure G\d = (<Tf +a±)/L. The account 
for rotation of the magnetization in the DW leads to a negative correction to this conductance. Our 
results differ from the results in papers published earlier. 



I. INTRODUCTION 



In ferromagnetic metals not only the charge of the electron but also the spin plays an important role in transport 
phenomena. A famous example is the observation of the giant magnetoresistance in magnetic multilayers, which can 
be explained in terms of a spin dependent electronic scattering. 

The presence of a domain wall (DW) in a ferromagnet can also change transport properties and this has been 
observed in a number of experiments. At the first glance, experimental data seems to contradict to each other. In 
Refs. |l],|2|,|],[§] it was found that the resistance of ferromagnetic wires and films decreases when increasing the external 
magnetic field, whereas in Refs. the resistance at zero magnetic field was found to be smaller than the one 
measured at high magnetic fields. 

In order to give a quantitative description of these experiments not only the DW contribution to the magnetore- 
sistance (MR) should be taken into account but also other mechanisms, as the anisotropic magnetoresistance, which 
arises due to the spin-orbit scattering size effects and the Lorentz contribution inside the domains. The main 

experimental difficulty in determining the DW contribution is to exclude the other effects. For example in Rcf. Q| 
the negative MR observed in Co films was interpreted in terms of DW scattering. However in Ref. || it was claimed 
that the predominant contributions to the observed magnetoresistance of Co films can be explained by a specific 
micromagnetic structure, which consists of stripe-domains with magnetization out-of-the-film plane. In addition, the 
films show closure caps at the surfaces with magnetization in plane and parallel to the current. Thus, the resistivity 
anisotropy might play a fundamental role. 

Understanding the details of these experiments is an interesting task. However, before taking into account all 
material specific characteristics of the experiments one should be able to describe general properties of electron 
scattering on domain walls. In this paper, we do not try to give an explanation of all these experiments, but solve 
an idealized model that may capture the most general features of transport in the presence of a DW. We calculate 
the resistance of a ferromagnetic wire with a DW and restrict ourself to the case when the magnetization of the 
ferromagnetic structure remains always perpendicular to the current. This assumption simplifies the situation because 
in this case anisotropic effects do not contribute to the change of the resistance. 

The DW contribution to the conductance has been considered in several theoretical works, in which different 
approaches have been used. For example in Refs. JlO|,[yJ quantum effects (weak localization) were taken into account. 
It was shown that a DW contributes to the decoherence of electrons leading to a decrease of the resistance. These 
effects may be important at very low temperatures when localization effects start playing a noticeable role. 

At higher temperatures the weak localization effects are not important and one may try to describe the magne- 
toresistance in terms of classical motion. In recent works, Refs. p^fDl, an increase of the resistance due to a DW 
was predicted on the basis of a Boltzmann equation. However, the collision term describing scattering of conduction 
electrons on impurities was introduced phenomenologically. 

The classical DW resistance was calculated also in Ref. [jL4| . In that work, it was shown that the DW resistance 
could be both negative and positive depending on the difference between the momentum relaxation times r-j- j for the 
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different spin directions. However, the cfassical Drude expression for the resistivity was used, in which the relaxation 
times r-[ i were introduced again as phenomenological parameters. In Refs. fH]Jl(| the resistance of a DW located in 
a point contact was calculated. 

The purpose of this paper is to derive a proper kinetic equation for the distribution function from a microscopic 
model and to calculate the DW resistance on its basis. We employ a standard approach based on microscopic equations 
for the quasiclassical Green's functions in the Keldysh technique. Assuming that the impurity scattering potential u s 
is spin dependent, we derive the kinetic equation for the matrix (in the Nambu and spin space) distribution function. 
As a result we come to the kinetic equation for the ditribution function / that is a 2 x 2 matrix in the spin space. 
The impurity scattering potential which enters the collision integral is also a matrix and this makes the equation 
considerably more complicated than the standard one that could be written for a spin independent scattering. 

Throughout this article we assume that the magnetization remains always perpendicular to the current. First we 
solve the derived kinetic equation in two simplest cases: a single domain and two domain structure with an abrupt 
DW (i.e. the width of the DW equals zero). In the case of a finite DW width we solve the kinetic equation assuming 
the potentials itf and u± do not differ much from each other. Even in this limit, it is hard to obtain analytical 
formulae for an arbitrary width of the DW. Two different limiting cases naturally arise and this allows us to obtain 
a solution for the distribution function. The first limit corresponds to a sharp DW (to a small exchange energy J ). 
The second limit corresponds to a smooth DW ( to a large J ). We note that only the second limit was analyzed in 
Refs. |lj,|l3|,[l7| . As in Refs. |lj,|l3|,[l7| we obtain that the DW increases the resistance of the system. However our 
formulae for the contribution of the DW to the resistance differ essentially from those presented in Refs. p^JT^ , [l7|| . 

The paper is organized in the following way. In the next section we introduce the model and derive the kinetic 
equation for the distribution function in a ferromagnetic wire neglecting quantum effects. We start from the mi- 
croscopic Hamiltonian (^|) with different scattering rates at impurities for spin-up and spin down- electrons. In the 
subsequent sections we calculate the conductance of the system in the diffusive limit. In section [II A we consider 
the case of a sharp DW when J <C D/w 2 , where D is the d iffusion coefficient, w is the width of the DW and J is the 
exchange field acting on the electron spin. In section [II B we calculate the conductance of a "slowly" varying DW, i. 
e. we consider the case D/w 2 <C J ■ It turns out that in the first case the conductance is always smaller than in the 
adiabatic case. In the last section we summarize our results. 



II. KINETIC EQUATION 

In this section we derive the kinetic equation for the matrix distribution function / starting from equations for the 
quasiclassical Green functions. The function / is a 2 x 2 matrix in the spin space. We assume that the impurity 
scattering rate depends on the spin directions but, for simplicity, we neglect such spin- flip processes as the spin-orbit 
interaction or the scattering by magnetic impurities. So, in our model each impurity scattering vertex is a matrix 
that does not commute with / and therefore the elastic collision integral has a nontrivial form. This fact has been 



ignored in Refs. [|12 17 1, where the collision integral was written phenomenologically. 

Using the derived kinetic equation we calculate the conductance of a mesoscopic structure which consists of two 
reservoirs and a ferromagnetic wire (or film) connecting the reservoirs (see Fig. 1). A domain wall is assumed to 
be present in the ferromagnet. We consider the diffusive limit, which means that the mean free path is the shortest 
length (apart from the Fermi wave length) in the problem. We solve the kinetic equation assuming the smallness of 
the parameter /3, defined as 

(1) 

where tr^j are conductivities for different spin directions. A precise relation between the conductivities cr-^j,, as well 
as the diffusion coefficients -Dj.,j, the corresponding scattering rates will become clear below. 

The assumption (3 <C 1 is valid for ferromagnets with the exchange energy J much smaller than the Fermi energy. 
We will consider two limiting cases: a) J << D^i/w 2 and b) J >> D^i/w 2 , where D^i is the diffusion coefficient 
for electrons with up and down spins, w is the width of the DW. The case a) corresponds to a sharp DW. The 
conductance in this case is smaller than the conductance of the structure without the domain wall. A finite width 
of the domain wall leads to a posit ive c orrection to the conductance. The second case corresponds to a smooth 
(compared to the magnetic length \J~D~Jj ) DW. In the limit of a large w the conductance of the structure is close to 
that of a structure without a DW. With decreasing the width of the DW, the conductance of the structure decreases. 
Our results significantly differ significantly from the results obtained in other works, where either the collision term 
was oversimplified Jl2|,0 , or the kinetic equation was not treated in a correct way |l3| . 

We choose the Hamiltonian of the ferromagnet in a simple standard form 
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H = J2 J dr{^{r)[-V 2 /2m + eV(r)- Jn- a]<tp s ,(r)} + H imp . (2) 

s,s' 

where V(r) is a smoothly varying (over the wave length \p) electric potential, J is the exchange energy, n is the 
unit vector directed along the magnetization orientation. The term Hi mp (r) describes the interaction of electrons 
with impurities and we assume that it depends on the spin direction. The origin of this dependence can be either 
the band structure or the intrinsec spin dependence of the impurity scattering potential [ fL2[ . If the magnetization is 
aligned along the z-axis, this interaction can be written as 

i 

As in Ref. |Q, we introduce new operators 

( ip s ,n = 1 , . 

^•• = \^+,n = 2 ' (4) 

In terms of the operators ipn tS and in the case of an arbitrary angle a between the magnetization vector and the 
z-axis the Hamiltonian (^) can be written as 



Himp = J2i, n , s J dr "^t,s( r ^ 3 ® &oU+(r - n) + (f ® a 3 cosa + f 3 ® a 2 sin a)u_(r - r,j).}^ n ^{r)} 



(5) 



where u± = (itf ± u\)/2, A = u-/u + and the matrix n is defined as h = T3 ® CT3 exp[— iaf3 ® o"i]. 

Introducing the operators VVs, Eq. (|4|), leads to an increase of the size of matrix Green functions written below. 
One has to deal not only with spin space but also with the Nambu one. Actually, this is not necessary if one 
consideres non-superconducting metals only. However, this extension of the size would become important if the metal 
wire we consider were in contact with a superconductor. Although we do not consider any superconductivity in the 
present work, we keep at the moment the Nambu space explicitly having in mind a possible generalization for the 
superconductivity. 

Now we define the Green functions in the Keldysh technique 

Gnn'iUA) = (1/0 < T C ^nAU)^At'k)) > ■ (6) 

where Tq means the time ordering along the Keldysh contour C. In a standard way we define the retarded (advanced) 
qR(a) anc j Keldysh Green function G as well as a matrix G composed of the matrices G R ^ and G (see e. g. Ref. 
|^9| ). One can obtain an equation for the matrix G in the usual way by summing the ladder diagrams in the cross 
technique (we neglect all crossed diagrams). This equation has the form 

{id t -H-S imp )G = l (7) 

where H = [— (l/2m)V 2 + V(t)]r 3 ® do — Jf 3 n, and the self-energy term Ei mp is given by 

S imp = n lmp u < G > u . (8) 

Here u = u+(l + An), < G >= v J d^ p J dfl/An ■ G and rii mp is the concentration of impurities. In the quasiclassical 
approach the density of states v is written in the main approximation with respect to the parameter J/ep, where ep 
is the Fermi energy. In this case v is the same for both spin-up and spin-down electrons. Notice that the r.h.s of Eq. 
(H) is a product of matrices, which in a general case do not commute. 
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FIG. 1. Upper Figure: The geometry considered in this article. The domain wall (DW) is situated in the middle of the 
ferromagnetic wire. Lower Figure: The z-component (n z (x)) of the magnetization in the F-wire. The solid line corresponds to 
a Bloch-like wall as calculated by Landau and Lifshitz. The dashed line corresponds to a linear DW 

In order to obtain an equation for the quasiclassical Green functions, we follow the standard way (see for example 
fl9f ): we write the equation conjugate to Eq. (7), multiply both equations by f 3 and subtract from each other. Then, 
we integrate the final equation over the variable £ p = vf(p — Pf) and obtain 



f 3 d t g + d t gf 3 + i(eV(t)g - geV(t)) + (v F V)g + i,J[n,g] = -(l/2r)(m < g > mg - gm < g > fh) 
We have introduced the quasiclassical Green function in the usual way 



g = (i/n)f 3 / d£ p G 



(9) 



(10) 



The matrix fh is equal to fh = 1 
of the matrix g are g R ( A '> and g 



Xn and r 



Vtli 



\u + \ 2 is the mean momentum relaxation rate. The elements 



g 



g 



~A 



(11) 



Eq.(|9|) is valid in a rather general case. In particular, it can be employed in the case of a superconductor- ferromagnet 
structure when the superconducting condensate penetrates into the ferromagnet. We use Eq.(|^) for a normal case, 
i.e. for F/S structures when one can neglect the penetration of the condensate into the ferromagnet F or for F/F' 
structures. In order to obtain the kinetic equation for the distribution function in the normal case, we represent the 
Keldysh component in the usual form 119| 



(12) 



where g R ( A > = ±t 3 <g) &o and / is a 4x4 matrix, whose elements are the components of the distribution function in 
the Nambu and spin space. This matrix can be represented in the form 



/ = f TQ + /3T3 



(13) 



The components f and / 3 are matrices in spin space. In the absence of spin-dependent interactions they are 
diagonal and related to the distribution functions for electrons n and holes p as follows: f = (1 — (n s + pj)) ■ 00; 
f 3 = ~(n s — pj) ■ (To. Taking into account Eqs.([H)|jL|), one can easily get from Eq.(||) the kinetic equation for the 
matrix distribution function 
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V)/ + Ut 3 [ft, /] = -(l/2r) [to 2 / + /to 2 - 2to < / > to] 
According to all previous definitions we can write 



(14) 



(1±A)2 



(15) 



and hence define o> ^ and -Df i without using any phenomenological approach. In our model the conductivities <7f i 
are equal to: er^ = e 2 vD^i = e 2 vD{l ± A) -2 , where D — v 2 t/3. Note that in the absence of superconductivity 
the distribution function has diagonal in the Nambu- space, and therefore one can take the component (1,1) of Eq. 
( |l4| ) and obtain 



(«fV)/ + iJ 



-(l/2r) 



TO 2 / + /TO 2 



2to < / > rh 



(16) 



where all matrices are now 2x2 matrices in the spin space. In particular, to — 1 + An and ri = (T3 exp[— ia&i]. Note 
that the left-hand side of Eq. (16) coincides with the left-hand side of the well known kinetic equation derived for a 
magnetic material (see for example Ref. pp[ | , where the kinetic equation is presented for a dynamic case in the absence 
of scattering by impurities). The solution / of this equation coincides with the component (1,1) of the distribution 
function /, which satisfies Eq. (|l4|). Since in this article normal materials (no superconductors) are considered, we 
will analyze Eq. (pf). One can exclude the spatial dependence of the matrices n and m performing an unitary 
transformation defined by 



/ = U ■ f ■ U + , U = (ToCOsa/2 + ia\ sina/2 . 



In this case one obtains an equation for the distribution function / 



(v F V)f + i{v F /2)a'{x) 



01,/ 



iJ 



= -(1/r) 



/- < / > +A 



5-3,/- < / > 



+ A 2 / - 0-3 < / > 03 



(17) 



The left-hand side of this equation differs from the one derived in Ref. jl3|]. In the latter there is an additional 
term of the form (a' (x)/Arri) [&i , {dF / dx)] + which, as we have shown, does not appear in the quasiclassical approach. 
Moreover, due to this term the kinetic equation of Ref. |l3|] violates the particle number conservation and therefore 
leads to wrong results. Notice also, that the collision term ( right-hand side of Eq. ( |l7| ) ) after the unitary rotation 
may not be diagonal in spin space. This fact was ignored in Refs. jl2];[L7]]. We will see in the next sections that in the 
case (D/w 2 ) <§; J, it is convenient to work with Eq. (|l7|) , while in the opposite case it is easier to solve the kinetic 
equation in its original form Eq. (|l6|). We assume that the system is diffusive (this implies the condition Jt << 1) 
. In this case one can expand the distribution function / in spherical harmonics and consider only the first two of 
them 



(18) 



where \i = cos9 and 9 is the angle between v F and the x -axis. Using Eqs. (16) and (|18|), one obtains two equations 
for the functions s and a 



vf9 x s + iJ [n,a] = — (1/2t) (m 2 a + am 2 ) 
(vp/3)d x a + iJ [n, s] = — (l/2r) (fh 2 s + sto 2 — 2tosto) 



(19) 
(20) 



In the second equation we have performed an averaging over the angle /1 . The boundary conditions at the interfaces 
with the reservoirs are given by imposing the continuity of the symmetric part s(x) of the distribution function (we 
assume a perfect contact of the F wire with the reservoirs) 



s(L) = tanh — do 



and 



s(0) = tanh ■ 



eV 



'IT 



-<7 



(21) 



(22) 
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Once we determine the distribution function /, we can calculate the current density using the following expression 

j = -jev^Y J deTra . (23) 

In the next sections we determine the resistance of a domain wall with a finite width. Here on the basis of 
Eq. (|l6|) , the conductance of a F'/F/F' mesoscopic system is calculated in the simplest cases: a single domain in the 
ferromagnetic wire and a two-domain structure in the F wire with an abrupt domain wall (i.e., w = 0, see Fig.l). In 
this case (the magnetization is parallel or antiparallel to the z-axis) , both parts of the distribution function s and a 
are proportional to <jq 3 . Therefore the commutator on the left hand side is equal to zero. FromEq. @ we find 

a = —m~ 2 vpd x s . (24) 

We substitute this expression into equation (poj). Taking into account that the right-hand side is zero, we obtain 
after integration 

s = s(0) + m 2 Ix/D , (25) 

The integration constant or, in other words, the "partial current" per unit energy I is found from the boundary 
condition ( p2| ) 

i=-jm- 2 F_a , (26) 

where m~ 2 = [1 + A 2 — 2An]/(l — A 2 ) and _F_ = tanh — tanh Substituting this expression into Eq. (p3|), 

we find the current and the differential conductance G = dl /dV\ v=0 

Gid = (2a/L)(l + A 2 )/(l - A 2 ) 2 = G T + G x (27) 

Here G-j- j = a^jL, and a — e 2 vD. Thus, the conductance has the usual form. We note that in terms of A the 
conductivities a^i are given by <7j,j_ = c/(l ± A) 2 , and hence the coefficient (3 defined in Eq. (|l|) is related to A via 
the relation f3 = -2A/(1 + A 2 ). 

Let us consider the same system with two domains in the F wire and with an abrupt DW located in the middle of 
the wire. In this case a = in the interval < x < L/2 and a = ir in the interval L/2 < x < L. Eqs. (|T9|-p0|) are 
solved in the same way as for the single domain case. For the symmetric part of the distribution function we obtain 

„ M f ootanh^tf^ +m 2 (Q)ix/D 0<x<L/2 , , 

S[X) ~ { s(L/2) + m 2 (ir)i(x-L/2)/D L/2 < x < L ' [ ' 



where rh (0) = m | a= o . The integration constant again is found from the boundary condition (P2|). We get for I 

I = a DF-/(l + \ 2 )L (29) 

and for the conductance 

G 2d = (2a/L)/(l + A 2 ) = 4G T G X /(G T + G L ) (30) 

This result has been obtained earlier (see Ref. (2) and references therein). In the next section we calculate G for the 
case when the magnetization (or the vector n) rotates in the y-z plane over a finite length w. 

III. CONDUCTANCE OF A DOMAIN WALL 

The problem of calculating the conductance for a system with a finite width of a DW is rather complicated. In 
order to simplify it, we make an assumption that the scattering times r^j, are close to each other, i.e. 

A<1. (31) 

This condition is met in fcrromagnets with an exchange energy J smaller than the Fermi energy. We consider again 
the system shown in Fig. [l]. The total length of the ferromagnetic wire is L. A Bloch-like DW is situated in the region 
(L — w)/2 < x < (L + w)/2 and separates two domains with opposite magnetizations. Thus, the effective width of the 
DW is w. It is not easy to obtain the exact solution of Eqs.(19|,|o|). However one can assume that condition ( |3i"| ) is 
satisfied and expand the functions s and a up to terms proportional to A 2 . We distinguish two cases: a) J <§C D/w 2 , 
which corresponds to a sharp DW; and b) J ^> D /w 2 . 



G 



A. Small exchange energy 



If the exchange field is weak (J <C D/w 2 ) or the DW wall is very sharp, one can easily solve Eqs. (fi9| 20 ) for 
an arbitrary form of the DW. We assume that the DW width exceeds the mean free path but is smaller than the 
magnetic length £j = y ' D / J . In this case, we expand the solution of Eqs. (|l9|,^0|) in the small parameters Jw 2 /D 
and A . In the zero order approximation, we get 



-ld x s 



and 



e + eV 

Dd x s = Iq] s = <To tanh — — h I x/D , 

where the 'partial current' is found from the boundary condition (^2|) and is equal to 

i = -(tqDF-/l 

In the first approximation we find from Eq. ( |l9| ) 

<ii = —ld x si ~ 2Xh(x)do 
The solution of Eq. ( po|) for the symmetric part Si has the form 

si = Iix/D + 2\i [ dxih{x x )/D 
Jo 

This and next corrections should satisfy zero boundary conditions. Therefore we find for I\ 



h 



-2XI <n> L 



(32) 

(33) 

(34) 
(35) 

(36) 
(37) 



where < (...) >l= 1/L J q {...)dx. As it follows from Eq.(|2|), the first correction does not contribute to the current. 
The zero order correction leads to the expression for the conductance given by Eq. (|3(]) if we expand it in the small 
parameter A (the case of an abrupt DW). In order to find a correction to the conductance due to a finite width of 
the DW, one has to find the second order corrections. One can see from Eq.(|23|) that only components of a,2 or §2 
proportional to &o contribute to the current. Therefore we take the trace in the spin space from Eq.(|l9|) and Eq.pOj) 
and find easily 



Tra 2 = -(l/D)TrI 2 



and 



where 



Trs 2 = -(2A) 2 Tr/ < n > L 



dxifi(xi)— < n > x 



ID, 



(38) 



(39) 



Tr/ 2 = -(l/2)Tr^ A 2 [l-4<n>|] . (40) 
Using Eqs. (CM) and Eqs.(f40|) we obtain the expression for the conductance which can be represented in the form 



G = G 



2,1 



1 + (2A) 2 -Tr <h> 2 L 



(41) 



This formula determines the conductance of the system under consideration for the case when the precession frequency 
J is smaller than the inverse time of diffusion of an electron through the DW. One can see that in the case of a 
DW with a finite width the conductance is larger than in the case of a sharp DW (c/. Eq. ((3^)), but smaller than 
the conductance in the single domain case (c/. Eq. (pTj)) . Note that Eq. ( f4l|) has been obtained in the limit of 
small A. Therefore the conductance Gid should be expanded in A (see Eq. (|30D) and terms of order higher than 
A 2 should be neglected. There is an interesting consequence from the result of Eq. (pl|). Let us consider the 
case of two DWs separating three regions of length d with homogeneous magnetization. For simplicity we assume 
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that the shapes of the DWs are described by a piece-wise linear function, which is characterized by a wave vector 
Q = (w/ir, 0, 0). If one defines the chirality vector as Sv c h = n(x) x n(x + 5x), where n(x) is the unit vector directed 
along the local magnetization, two cases should be distinguished: a) the DWs have different chirality. In this case 
Tr < n >\= (2/L 2 )(d 2 + I6(w 2 /ir 2 )). Thus we see that an additional DW will decreases the conductance of the 
system, b) the chiriality vectors have different signs. In this case Tr < ft >\= (2/L 2 )d 2 , and hence the contributions 
of both DWs to the conductance cancel each other. This result can be generalized easily for an arbitrary number of 
DWs 

Now we calculate the spatial distribution of the electric field in the ferromagnetic wire showed in Fig. [j]. The 
electric potential V(x) is given by the expression (see, for example, Ref. p|) 

V(x) = (l/4)Tm / des (42) 

According to Eqs. ( |33| ) and ( |39| ) the electric field E(x) = —d x V(x) in the ferromagnetic wire is given by 

E(x) = (V/L){1 + (2A) 2 [< cosa > L (< cosa > L - cosa)+ < sina > L (< sin a > L - sin a)]} . (43) 
For example, if we consider the structure of the Bloch wall which has been calculated by Landau and Lifshitz [ p2[ 

cosa = tanh[(a; - L/2)/w]; s'ma = cosh _1 [(a; - L/2)/w] (44) 

we obtain 

E(x) -E = E {2\) 2 {itw/L) {(ttw/L) - cosh" 1 !^ - L/2)/w}} , (45) 

where Eq = (V/L). In Fig.^ we plot the dependence E(x) given by Eq. ([45|). One can see that in the region of the 
DW the electric field decreases; this means an increase in the local conductivity. 

In the next section we consider the case of a strong exchange field or of a wide wall, i.e. the case w ^> £j . 




B. Large Exchange Energy 

Now we consider the case with a large exchange energy or a slow variation of the direction of the magnetization 
within the DW (w » y/D/J). In this case, the problem becomes more complicated because we cannot neglect the 
commutator on the left hand side in Eqs. (|Tj^,pO|) and cannot find a solution of these equations even for the case of 
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small A. Therefore we simplify the problem assuming that the shape of the DW is described by a piece-wise linear 
function (see Fig. [|) 



03 in region I 

63 exp{-i(Tr/w)(x - (L - u))/2)(Ti] within the DW (46) 
— a 3 in region III. 

Obviously the results for other shapes of the DW like that given by Eq. (Q) will differ from ours only by a numerical 
factor. We again expand the solution in the small parameter A, i.e. a = &q + d\ + 0,2 + .... The zeroth order terms 
can be obtained easily as before and they are given by Eqs.([32],[53|). The first correction a\ is given again by Eq. ( |35| ) 
and the first correction for the symmetric part s\ obeys the equation 

Dd 2 xx h - U[h{x), si] = 2XI d x h . (47) 

This equation can be solved for the case of h(x) given by Eq.(46|) with the help of an unitary transformation (a 
rotation in the spin space). We do not need to find the second order corrections hi and S2, since the sought correction 
to the conductance can be expressed in terms of s\. Indeed, let us write the equation for Trs2 which follows from 
Eq.© 

Tr{Dd x s 2 + (v F /3)[X 2 a + 2\ha{\ - J 2 } = , (48) 

where I2 is the integration constant which is related to 0,2 : Trf7oa2 = —II2/D. We integrate this equation from to 
L taking into account the boundary conditions at x = and x — L: §2 = 0. After simple transformations we obtain 

Tr{/ 2 - 3\ 2 I }/D = 2ATr<7 { / dxs^nix)}; . (49) 

Jo 

As before, the correction to the conductance is determined by 1 2. We see that in order to find this correction, one has 
to solve Eq. ( ^7| ) for §\. This equation can be solved with the help of the unitary transformation 

h = U ■ S ■ U + , U = <T cosa/2 + i<Tisina/2 (50) 

This rotation transforms the vector n into N = 03. Performing the [/-transformation, we obtain instead of Eq.(|49|) 

Tr(7o{/2 - 3X 2 I }/D = 2Xd x aTr & a 2 < S(x) > . (51) 

After the [/-transformation Eq.(p7|) acquires the form (in the region of the DW) 

d 2 xx S - (Q 2 /2){S - axSax) + iQ{atdJ - d x Sa x ) - iJ[N, S] = 2Xa 2 I Q/D , (52) 

where Q = d x a and N = 03. This equation is valid in the region of the DW, whereas in the regions I and III we 
have to set Q = and to take into account that in the region III N = —as. The matrix S should be represented as 
a sum: S = S\&\ + S f 2 (5'2 + The components Sk are given by a linear combination of the eigen-functions of Eq. 

([52|). They obey zero boundary conditions at x — and x = L and should be matched at x = L\ = [L — w)/2 and 
x = L 2 = (L + w)/2. The eigen-values of Eq.^) (Sk ~ exp(ra)) are determined by the equation 

K 2 {K 2 + Q 2 ) 2 = {K 2 ~Q 2 )iej , (53) 

where ^J 2 = 2J/D. In a general case a solution of Eq. has a cumbersome form. We represent here the form of a 
solution for Tro^^x) in the region of the DW which we are interested in: 

Tr& 2 S{x) S -2XTrJ (£jQ/D) Im[exp(-(1 + i)(x - L x )Hj <&) + exp((l + i)(x - L 2 )/ti V^)] . (54) 

We dropped terms of the higher order in the parameter Q£j ~ £,j/w. Using this expression and Eq. (pl|), we readily 
get the expression for the current and for the conductance G 

G=Gid { i ~i3 x2 ) ' (55) 
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where G\d is the conductance for a homogenous magnetized wire (see Eq. (p7|)). Again terms of order higher than 
A 2 should be neglected. Note that G\d is always larger than the conductance in the case of a two domain wire Gid 
(see Eq. (|30|)). Eq. ( |55| ) shows that the DW decreases the conductance compared to the conductance Gid of a single 
domain F wire. Our result is sketched in Fig. ||. We see that within our approach a DW with a finite width is always 
a source of resistance. 



G ld 



G 2d 

(in arbitrary units) 

FIG. 3. Schematical representation of the conductance as a function of the width w of the DW. In the intermediate region 
(dashed line) the curve is extrapolated from our results. 

IV. CONCLUSION 

Using a simple microscopic model (equal density-of-states but different impurity scattering times r^j. for electrons 
with spin up and down), we have derived the kinetic equation for the matrix distribution function. The derivation has 
been performed by a standard method on the basis of microscopic equations for the quasiclassical Green functions in 
the Keldysh technique. This equation can be applied to the studies of transport in, for example, ferromagnets with 
a non-homogeneous magnetization. 

We have employed this equation to calculate the conductance G in a mesoscopic F'/F/F' structure. We have 
assumed that the parameter A = (t^ — Tf)/ 2(r^ + T|) is small and the length of the F wire L is shorter than the spin 
energy relaxation length. Two different limits appear which are determined by the product of the exchange energy 
J and the diffusion time t w = w 2 /D of electrons through the DW. In the limit t w J <C 1 and a very thin DW the 
conductance of the structure (per the unit cross-section area) is equal to Gid = 4<7-r<7j(<7j + a\)/L. The account for 
a finite width of the DW leads to an increase in the conductance by a normalized amount of the order (Xw/L) 2 . We 
have also calculated in this limit the spatial distribution of the electric field in the F wire. The electric field has a 
minimum in the center of the DW which corresponds to an enhanced local conductivity. In the other limit t w J » 1 
(adiabatic variation of the magnetization in the DW) the conductance coincides in the main approximation with that 
of a single domain structure G\d = (c| + <7jJ/L. The account for rotation of the magnetization in the DW leads to a 
negative correction to the conductance of the order -\ 2 (t w J)~ 3 / 2 (w/L). Our results differ from those published earlier 
p2| , ^3|Jl7| because in the latter works the collision term was written phenomenologically. In particular the matrix 
character of the impurity vertex was not taken into account. 

We would like to thank SFB 491 for financial support. 
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